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1. Introduction 



Fraunhofer diffraction by a circular lens or aperture is 
a ubiquitous phenomenon in optics in general and ra- 
diometry in particular. Figure 1 illustrates two practical 
situations in which Fraunhofer diffraction occurs. In the 
first example, diffraction limits the ability of a lens or 
other focusing optic to focus light. According to geo- 
metrical optics, it is possible to focus rays incident on a 
lens to converge at a focal point. In practice, even with 
the focal point at the center of a finite-sized circular 
detector, some of the light incident on the lens fails to 
be collected, giving rise to a "diffraction loss." 

The second example can occur when using a source, 
such as a blackbody cavity source viewed through a 
small "pinhole" aperture, to calibrate optical systems 
intended for observing celestial objects or for use in 
remote- sensing applications. In this case, the angle sub- 
tended by the cavity opening at the aperture is larger 
than the angle subtended by the detector pupil on the 
dark side of the aperture [1]. Because the former angle 
is larger than the latter one, geometrical optics suggests 
that the total power detected depends only on the black- 
body temperature and a geometrical factor related to the 
pinhole aperture, detector optics, and relative separa- 



tion, because the detector pupil is overfilled. However, 
diffraction leads to losses in the total power reaching the 
detector. 

All of the above diffraction losses have been a subject 
of considerable interest, and they have been considered 
by Blevin [2], Boivin [3], Steel, De, and Bell [4], and 
Shirley [5]. The formula for the relative diffraction loss 
in spectral power is already well known, and it is noted 
below. However, a formula for the diffraction loss in the 
total power for the case of a Planckian source such as a 
star or blackbody appears to have been given only in the 
high-temperature limit. This article reports a formula 
for the diffraction loss in the power for such a source at 
all temperatures. The formula is useful for predicting 
the loss in cases such as the examples discussed, and it 
can be used as an independent check for numerically 
calculated diffraction losses. 

For the radiation present at a given wavelength A , the 
diffraction loss can be described in terms of a unit-less 
parameter, v = Iti^RIK. R is the radius of the lens or 
aperture. The angle ^''is defined in either of two ways: 
either 2^ is the full angle subtended by the detection 
pupil at the focusing optic, or 2^ is the full angle 
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Fig. 1. Two optical systems are shown, for which the (spectral) power detected is subject to 
Fraunhofer diffraction. In the first system (a), a lens concentrates radiation on the detector. In the 
second system (b), radiation from a blackbody cavity may pass through the blackbody cavity 
opening, pass through the pinhole aperture, and reach the detector pupil. 



subtended by the blackbody cavity opening at the 
pinhole aperture. For a source in thermal equilibrium 
at temperature T, the effects of temperature on the 
distribution of spectral power enter relevant equations 
through the ratio, C2l{X.T)=Av. Here, C2 = hclk = 
1.438 7752(25) X 10"^ m K is the second radiation con- 
stant, where the number in parenthesis is the one- stan- 
dard-deviation uncertainty in the last two digits. This 
implicitly defines another unit-less parameter, 
A = CjlilT^'^RT). The relative effects of diffraction on 
the total power only depend on this one parameter, A . 
As a final note, if the detection pupil subtends a 
considerable angle, the diffraction losses can be approx- 
imately given as the weighted average of losses arising, 
in the case of an infinitesimal detector, for values of ^ 
varying between ^ minus the half-angle subtended by 
the detector and ""f^plus the same half-angle. Analogous 
weighting can be used, in the context of the first exam- 
ple, to account for the finite angular diameter of a source 
(smaller than 2^). Such averaging is discussed else- 
where [5] but is not discussed further in this work, 
which henceforth assumes a single effective value of ^. 



2. Derivation of Formula 

In the Fraunhofer case, the ratio of the true spectral 
power to the spectral power predicted by geometrical 
optics is 



F(v)=l - J^(v) - J?(v). 



(1) 



Here, i„iv) is a cylindrical Bessel function. It is neces- 
sary to incorporate this wavelength-dependent factor 
when evaluating the ratio of the total power detected to 
the total power predicted by geometrical optics. For a 
source whose spectral output obeys the Planck distribu- 
tion, this unit-less ratio is given by 



F{A): 



;dvv^[exp(Av) - 1]-'[1 - J?(v) - J?(v)] 



;dvv'[exp(Av) - 1]"' 



(2) 



The value of the denominator is familiar, being given 



by 
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/dvv''[exp(Av) — 1] ' = 2/dvv^^ exp(— nAv) 

n=10 



(3) 



Here, ^(« ) is the Riemann zeta function. When evaluat- 
ing the numerator, two techniques have been found to be 
helpful: one for the "low-temperature" case (defined as 
A > 4) and one for the "high-temperature" case (defined 
as A < 4). By such definitions, "high" and "low" temper- 
ature cases arise depending on the average diffraction 
loss for a source at temperature T. This depends on T 
and the geometry through A . In either case, one first 
makes use of the relation. 



Jo[2v sin(0/2)] = J?,(v) -I- 2 2 Ji(v)cos(m0), 



from which follows. 



(4) 



J?(v) + Jf(v) = — Jddil + cos0)Jo[2vsin(0/2)]. (5) 

iTT 



Low-Temperature Case 

To evaluate the numerator in the low-temperature 
case, series expansion of the latter Bessel function in Eq. 
(5) yields 

J^(v)H-J?(v) = - i '■"/''''' ' /d0[l - sin'(0/2)]sin'X0/2) 



TT s^{) {S !)■ 



(-I)'..-' 



^^,5 ,,!)■ 



(2.y-l)!! (2iH-l)!! 



= 22 



(-l)'v^' 



.=0 (.v!)^ 



. (2.V)!! 
(2^- 1)!! 



(2s + 2)\\ 



{2s + 2)\\ 



(6) 



This permits one to rewrite Eq. (2) as follows: 
e(A)^6A-t(4)[l -F(A)] 

= 2 I dvv'exp(-«Av)[J^(v) + J?(v)] 



= 2l 



= -s 



{-\)\2s + 'i)\^{2s + A) 



{s\)W 



(2.9 - 1)!! 
(2^ + 2)!! 



{2'Tit*\2s + l,)\{2s - l)!!fi2„+4 
{s^:f{2s + A)\{2s + 2y.\A^'*^ ' 



(7) 



In this expression, 5, is a Bernoulli number. The sum 
converges for all A > 2, and one obtains an error in F(A ) 
estimated to be as small as 10 ''' for A > 4 if one sums 
up to 5 = 28." 



High-Temperature Case 

To evaluate the numerator in the high-temperature 
case, the relation in Eq. (2) is rewritten to render 

G(A) = 6A-t(4)[l-F(A)] 



with 



= 2;dvv'exp(-«Av)[J^(v) -i- J?(v)] 

n=10 



:2ei(nA), 



ei(A) = /dvv'exp(-Av)[J?(v) + J?(v)]. 



Use of Eq. (5) and the relation, 
/dvexp(-av)Jo(^v) = 



one has 

V dA/ 2-77 VA'-H4sin"(0/2) 
Making the abbreviation, w = 4/{A^ -n 4), one has 



V^^^' 



1 -H cost 



(8) 



(9) 



(10) 



(11) 



Gi(A) = 



'dA 



4 "''^ , , cos^rf) 

/d(/)- 



'dA 



ttV'A^h-4 Vl - wcosV- 

/ 4 [K(w) - E{w) 

1 IT Va'h- 



dVfV^ 



H-4 



dA/ [ 77 



{K{w) - £(w)] 



(12) 



where E{w) and ^(w) are respectively complete elliptic 
integrals of the first and second kind, defined according 
to the convention,' 



ir/2 




E{w)= jTd^Vl — vvsin^i^. 



Ttll 

Kiw)= f 



Vl — wsin^cj) 



(13) 



(In a different convention, w is replaced by w^ in the 
integrand but nowhere else.) From the properties of 
elliptic integrals [8], one may obtain 



Gi(A) = 



ttA' 



4(2-i-A')£'(w) - A^K(w) 



(14) 



and 
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Qi(A)- 



_d_ 
dA 



4Eiw) 



.TrAWA^ + 4. 



(15) 



The latter result is easily obtained, because Qi(A) has 
already been written as the third derivative of an expres- 
sion with respect to A . 

The closed-form result in Eq. (14) may be directly 
substituted in Eq. (8), after which one may sum over n . 
Because the summand can be integrated according to 
Eq. (15), however, the Euler-Maclaurin formula may be 
used to obtain a more easily evaluated, asymptotic ex- 
pression for 2(A): 



QiA) = ^+ 2 Q^A) + - ;dA'e,(A') + - QANA) 

TTA ,1=1 A NA Z 



'liki^'^-'iM^'^Q^^^i 



+ R„,iN,A), 
(16) 



where 5, is again a Bernoulli number. Here a new func- 
tion, 



e,(A) = G,(A) 



ttA' 



(17) 



has been introduced as the summand, so that the largest 
part of 2i(A ) is extracted from the sum prior to summa- 
tion. 

The parameters m and N are both positive integers. 
Rm{N,A), the remainder (i.e., error) term, usually de- 
creases initially with increasing m, but then rapidly 
diverges. It is therefore good to choose a suitable combi- 
nation of N and m , for a given value of A , to have an 
efficient yet accurate result. Increasing A^ should usually 
improve results, and one obtains an error in F{A) esti- 
mated to be as small as 10 ''' for A < 4 if one uses 
TV = 60 and m = 5. 

Explicitly, QiA) is given by 



4/(3) 'V-i. 
Q(A) = ^+ ^Q,{nA) 

TTA n=l 



2 j mN+l)+2a\N+2)]E(co)-a^K(co) 



1/1 1 



(A+aT 



^Tia'\ N'^3N' 3N''^5N'' 



1 -' B2S leAa)E((o)-a\ia)K((o)]f,, 
7,a' ii {2sy. N''-\4+aY'*'"' 



+ R,„(N,A). 



■N-l 



(18) 



One should include only the first m — 1 terms in the 
parenthetic expression preceding the second sum. Two 
abbreviations, a = NA and co = 4/(4 -i- a^), have been in- 
troduced. In the second sum, the first several constants 
are given by 



B2_J_ B4 

2! ~ 12' 4! ' 



1 B, 



1 



B, 



1 



720' 6! 30240' 8! 



1209600 ' 



/i = 4, /2 = 24, f, = 48, /4 = 2880. 



(19) 



Likewise, the first several e- and ^-functions are given 
by 

ei(a) = 96 -I- 68a" -H I9a\ 

ki(a)= \2 + la\ 

e2(a) = 5120 + 6112a" + 2880a" -1- eSOa" + 117a', 

k2(a) = 640 + 624a' + 216a* + 51 a\ 

e^a) = 1720320 + 2908160a" + 2093376a'* + 

8295520" + 196164a'' + 199170'° + 3726a'", 
k,(a) = 215040 -I- 316480a" + 193072a* + 628920*^ + 

8265a' -H 2046a'", 
e4(a) = 33030144 + 72310784a" + 70087936a* + 

394298880" + 14163680a' + 33501280'" + 

572853o'" -I- 21062o'* + 6121a'", 
k4(a) = 4128768 + 8135680a" + 6993408o* + 

3421440a" -1- 1037920o' + 2195160'" + 

8820o'"h-3601o'*. 



(20) 



3, Evaluation of Formula 

In the high-temperature limit, this gives a result of the 
form 



F(A)=1 



[4A^(3)/TT-K9(AMnA)] 
6^(4) 



(21) 



and evaluating the successive terms in Eq. (21) may be 
both difficult and unfruitful, because the series appears 
to be both very complicated and slow to converge. Note 
that the leading terms in Eq. (21) are consistent with 
other works cited, such as the work of Blevin. In the 
low-temperature limit, one has 



F(A) = 



IOtt" 77* 

21A" 4A*' 



(22) 



Figure 2 shows F(A) for a range of values of A 
bridging the low- and high-temperature regions. The 
solid curve indicates "exact" results obtained using Eq. 
(7) for A > 4 and Eq. (18) for A < 4. The dashed curves 
indicate approximate results for low- and high-tempera- 
ture limits, obtained using the respective formulas in 
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Ife. 




Fig. 2. F{A) vs A according to Eq. (7) and Eq. (18) (solid line), and 
Eqs. (21-22) (dotted line). The latter, approximate results apply in the 
limits of small A and large A . 

Eqs. (21-22). Likewise, Table 1 shows sample values of 
F(A) over a similar range. For < A < 0.2 and A > 8, 
inclusion of terms shown in Eqs. (21-22) yields an error 
in F{A) smaller than 0.001. While this immediate dis- 
cussion helps provide a sense of the behavior of F(A), 
accurate values should be found using Eq. (7) and Eq. 
(18). 
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Table 1. F{A) at sample values of A 



A 


F{A) 


A 


F{A) 


A 


F{A) 


A 


F{A) 


0.2 


0.9526 


2.2 


0.4643 


4.2 


0.2055 


6.2 


0.1077 


0.4 


0.9041 


2.4 


0.4259 


4.4 


0.1912 


6.4 


0.1018 


0.6 


0.8543 


2.6 


0.3907 


4.6 


0.1781 


6.6 


0.0964 


0.8 


0.8033 


2.8 


0.3587 


4.8 


0.1663 


6.8 


0.0914 


1.0 


0.7514 


3.0 


0.3297 


5.0 


0.1555 


7.0 


0.0867 


1.2 


0.6994 


3.2 


0.3034 


5.2 


0.1457 


7.2 


0.0824 


1.4 


0.6481 


3.4 


0.2797 


5.4 


0.1367 


7.4 


0.0784 


1.6 


0.5984 


3.6 


0.2582 


5.6 


0.1285 


7.6 


0.0746 


1.8 


0.5508 


3.8 


0.2389 


5.8 


0.1210 


7.8 


0.0712 


2.0 


0.5060 


4.0 


0.2214 


6.0 


0.1141 


8.0 


0.0679 
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